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I. INTRODUCTION 

The understanding of hadronic systems in terms 
of their elementary degrees of freedom have been, 
and still is, one of the most challenging problems in 
particle and nuclear physics over the last ten years. 
The phenomenological properties of hadrons are now 
rather well understood in terms of models, like the 
constituent quark model or the bag model. The un- 
derstanding of their properties from the original La- 
grangian of QCD is however still under active debate. 

In nuclear physics, the properties of nuclear struc- 
ture in terms of the exchanges of pions are also well 
known. They are described by using a phenomenolog- 
ical nucleon-nucleon potential expressed in terms of 
the exchanges of one pion, two correlated pions and 
so on. However, their complete description from an 
effective chiral Lagrangian is still missing. 

A common difficulty in both domains is the descrip- 
tion of relativistic bound systems. This description 
should be nonperturbative from the start in order to 
be able to find, for instance, the physical mass of the 
bound state from the pole of the scattering amplitude. 
The problem is especially acute when the interaction 
coupling constant is large. 

One of the most relevant approaches aimed at 
studying relativistic systems of interacting particles 
is light-front dynamics (LFD) proposed initially by 
Dirac [T]. LFD is a form of Hamiltonian dynamics 
which deals with the state vector defined not at a fixed 
time moment, but on the light front plane t-\-z = 0, in 



its traditional form. The state vector is then usually 
decomposed in a series of Fock sectors, each contain- 
ing a fixed number of particles. 

The use of LFD to investigate relativistic bound 
states has been advocated for a long time. However, 
while the dynamics of few-body systems, based on a 
phenomenologically constructed interaction, has de- 
veloped rapidly, application of LFD to field theory be- 
yond a perturbative framework is not yet under com- 
plete theoretical control. This is due to the fact that 
any practical calculation relies on taking into account 
only a restricted number of Fock sectors in the state 
vector decomposition or, in other words, on the Fock 
space truncation. This approximation strongly com- 
plicates the renormalization procedure, in contrast to 
that in standard perturbation theory. Indeed, the full 
cancellation of field-theoretical divergences which ap- 
pear in a given Fock sector requires taking into ac- 
count contributions from other sectors. If even a part 
of the latter is beyond our approximation, some di- 
vergences may leave uncancelled. Mathematically, it 
reflects itself in possible dependence of approximately 
calculated observables on the regularization parame- 
ters (e. g., cutoffs). This prevents to make any phys- 
ical predictions if we cannot control the renormaliza- 
tion procedure in one way or another. 

In a previous study [2] (see also references therein) 
we have developed an appropriate renormalization 
procedure - the so-called Fock sector dependent renor- 
malization (FSDR) scheme - in order to keep the can- 
cellation of field-theoretical divergences under perma- 
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nent control. Our approach is based on the covariant 
formulation of LFD (CLFD), where the state vector 
is defined on an arbitrary light-front plane character- 
ized by a light-like four-vector lo [3] and given by the 
equation cj-x — 0. The covariant formulation is nec- 
essary in order to control any violation of rotational 
invariance, including that which is caused by the Fock 
space truncation. In particular, this is important in 
order to formulate, in an unambiguous way, the renor- 
malization conditions one should impose on the bare 
couphng constant (BCC) to relate it to the physical 
one. 

In Ref. 12] we calculated the fermion state vec- 
tor and the electromagnetic form factors within the 
Yukawa model and QED in the lowest nontrivial ap- 
proximation, when the state vector includes only two 
Fock sectors given by one constituent fermion and 
one constituent fermion coupled to one boson. For 
this two-body Fock space truncation, the electromag- 
netic form factors are identical to those obtained in 
the second order of perturbation theory, giving rise to 
a Schwinger-type correction to the fermion magnetic 
moment. Note that this result is not surprising, in 
spite of the fact we have not done any expansion in 
powers of the coupling constant, since no other contri- 
butions to the fermion electromagnetic vertex, apart 
from the perturbative ones (resummcd to all orders) 
is generated in the two-body truncation. 

We shall present in this work the calculation of 
the fermion anomalous magnetic moment (AMM) 
within the same Yukawa model, but for the three-body 
Fock space truncation, when the state vector includes 
an additional Fock sector containing one constituent 
fermion coupled to two scalar bosons. The presence of 
three-body states gives rise to nontrivial nonperturba- 
tive contributions to the AMM, which can not be fully 
incorporated in perturbation theory. Besides that, the 
Yukawa model is a quite nontrivial one from the point 
of view of the renormalization procedure, since it ex- 
hibits simultaneously mass, vertex, and wave function 
renormalization. 

The plan of the article is the following. We recall in 
Sec. [n] the main features of our nonperturbative ap- 
proach and the renormalization procedure. The eigen- 
value equations are derived in the three-body trunca- 
tion in Sec. |III[ We calculate the fermion electromag- 



netic form factors in Sec. IV and present our numerical 
results for the AMM in Sec. |Vl Our conclusions are 
drawn in Sec. |VI| The appendices collect all necessary 
details to calculate the AMM. 



II. BOUND STATE SYSTEMS IN 
LIGHT-FRONT DYNAMICS 

A. General framework 

The state vector describing any relativistic sys- 
tem with total four-momentum p forms a representa- 
tion of the Poincare group. The four-momentum op- 
erator squared is one of the Casimir operators of 
this group and the state vector satisfies the equation 



(1) 



where M is the mass of the physical system under 
consideration and -p^ = M^. 

LFD serves as an effective and convenient tool to 
solve this eigenvalue equation. Indeed, one of the main 
advantages of LFD is that, due to kinematical con- 
straints, the vacuum state of a system of interacting 
particles coincides with the free vacuum, and all inter- 
mediate states result from fluctuations of the physical 
system. One can thus construct the state vector in 
terms of combinations of free fields, i.e. decompose it 
in a series of Fock sectors: 



E 



dDn(f>n{ki, . . . ,fc„;p) |n) 



(2) 



where \n) is the state containing n free particles with 
the four-momenta fci,...,fc„ and 0„'s are relativis- 
tic n-body wave functions, the so-called Fock compo- 
nents. The phase space volume element is represented 
schematically by dDn- In the following, we shall re- 
strict our study to a physical system composed of one 
fermion and n — I bosons. In that case 



\n) = a\h)c\k2)...c\K)\0), 



(3) 



where a' and are fermion and boson creation oper- 
ators, respectively, and 

0„(fci, . . . fc„;p) = u(fci)V'„(fci, . . . k„-p)u{p), (4) 

where u's are bispinors. To completely determine the 
state vector, we normalize it according to 



(5) 



With the decomposition ([2| , the normalization condi- 
tion ^ writes 



oo 

E 



In = 1, 



(6) 
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where /„ is the contribution of the n-body Fock sector 
to the full norm. Explicit formula for /„ can be found 
in Ref. 12]. 

In the following, we shall use CLFD |3j as a general 
framework. The covariance of our approach is due to 
the invariance of the light front plane equation. This 
implies that uj is not the same in any reference frame, 
but varies according to Lorentz transformations, like 
the coordinate x. It is not the case in the standard for- 
mulation of LED where w is fixed to uj — (1, 0, 0, —1) 
in any reference frame. 

The light-front momentum operator Pp can be con- 
structed from the energy-momentum tensor. It is de- 
composed according to 



Pn 



p(0) 

p 



P ' 



(7) 



where the two terms on the r.-h. s. are, respectively, 
the free (i. e. independent of the coupling constant 
and counterterms) and interaction parts of the four- 
momentum operator. The operator P™* is related to 
the interaction Hamiltonian on the light front 

by 



The vertex function F„ will be represented graphically 
by the diagram of Eig. [ij With the definition 



oo „ 

= 53 / dDnu{ki)Tn{ki,...,kn,p)u{p)\n) 



the eigenvalue equation ([T]) writes [2] 



(12) 



Qip) 



27r 



dr 



-h'-\ut) -g{p), (13) 

T 



where i^"** is the interaction Hamiltonian in momen- 
tum space: 



(14) 




(n — 1) bosons 



P 



Up / W\x)5{uj-x)d'^x. 



(8) 



Erom the general transformation properties of the 
light-front plane lo-x = 0, one can derive the following 
conservation law [3; for each Eock component: 



fci + fc2 



kn = P 



(9) 



where the quantity r„ is a measure of how far the n- 
body system is off the energy shell^. It is completely 
determined by the conservation law ^ and the on- 
mass shell condition for each individual particle mo- 
mentum: 



2uj-p r„ = (s„ - Af^), 



(10) 



where s„ = (fci + . . . + /s„)^. 

It is convenient to introduce, instead of the wave 
functions (j)m the vertex functions F„ (which we will 
also refer to as Eock components), defined by 



u(fci)F„u(p) = (s„ - 



2iVp Tn(j>n- (11) 



FIG. 1: Vertex function of order n for the Fock space 
truncation of order TV. 



With the form (T3|, the eigenvalue equation can 
thus be represented graphically, using the same rules 
as those derived in Ref. [3] for the calculation of matrix 
elements of the S'-matrix. This graph techniques was 
developed by Kadyshevsky |3j and transformed to the 
case of CLFD in Ref. [5]. 

The substitution of the decomposition ( |l2| ) into the 
eigenvalue equation ([Ts]) results in an (infinite) sys- 
tem of equations for the Eock components. In order 
to solve this system in practice, we should make it 



finite, i. e., truncate the decomposition (12), or equiv- 
alently ([2]), by retaining only those Fock sectors where 
the number of particles does not exceed some maximal 
value N . The finite system can be solved numerically 
and nonperturbatively, that is, for any value of the 
coupling constant. This approach was developed in a 
series of papers [H [BHH] • In Ref. [5] it was applied to 
the case of the two-body truncation, i. e. for N = 2. 



^ The term "off the energy shell" is borrowed from the equal- 
time dynamics where the spatial components of the four- 
momenta are always conserved, but the energies of interme- 
diate states are not equal to the incoming energy. 



B. Renormalization conditions 

In order to be able to make definite predictions for 
physical observables, one should also define a proper 



B Renormalization conditions 
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renormalization scheme which allows to express, like 
in perturbation theory, observables through the physi- 
cal coupling constant and masses and exclude the bare 
ones. The basis of the state vector decomposition, 
i. e. the states \n) in Eq. ([s]) is constructed from free 
physical fermion and boson fields, with their physical 
masses m and fj,, respectively. The interaction Hamil- 
tonian contains the corresponding mass counterterms 
(MC's) Sm and S^^ responsible for the fermion and 
boson mass renormalization. Since we will not con- 
sider here the fluctuations of the boson in terms of 
fermion-antifermion pairs, we have to set Sfj,^ — 0. 
The MC Sm is determined from the eigenvalue equa- 
tion ( 13 ) by demanding that the bound state mass M 



is equal, for the ground state, to the physical mass m 
of the constituent fermion. For this reason, we will 
distinguish M and m only when it is necessary. Oth- 
erwise, we will set M = m. 

Besides MC's, the interaction Hamiltonian includes 
also the BCC go. The latter is determined, as in per- 
turbation theory, by relating the on-energy-shell two- 
body vertex function r2 to the physical coupling con- 
stant g. As follows from Eq. (|9|, taking r2 on the 
energy shell is equivalent to setting T2 — 0. Once 
M is identified with m, the latter cond ition reduces 
to S2 = (fci + ^2)^ = m'^ [see Eq. (10)]. Below, for 
brevity, we will denote the on-energy-shell two-body 
vertex function as r2(s2 = m^) to indicate that its 
arguments are connected by the corresponding kine- 
matical constraint. 

Being a solution of the system of eigenvalue equa- 
tions (13), r2 depends on the BCC go. Hence, relat- 



ing r2 to g is equivalent to relating go to g, which just 
means coupling constant renormalization. This can be 
most easily done starting from the three point func- 
tion with all undressed on-mass-shell external lines, 
called r2. It is connected with the physical coupling 
constant by the following standard relation (see e.g. 
Ref. ID): 



Zb = .g. 



(15) 



where the Z-factors are the so-called field strength 
renormalization constants for the fermion (/) (both 
in the initial and final state) and boson (b) lines, re- 
spectively. This condition can also be recovered by 
demanding that the residue of the fermion-boson elas- 



full fermion self-energy E( 



(16) 



tic scattering amplitude at S2 



equals g . One 



can thus deduce the expression for Zf vn terms of the 



and similarly for Zf, as a function of the full boson 
self- energy. 

The two-body Fock component r2 being a solution 
of the eigenvalue equation ( 13 1 does not coincide with 
r2. By definition, r2 has no radiative corrections to 
any of its three legs, while r2, on the contrary, in- 
cludes such corrections. The relation between these 
two vertex functions taken off the energy shell is rather 
complicated in LFD. Fortunately, we need to know it 
on the energy shell only, where it simplifies strongly, 
because the on-shell radiative corrections mentioned 
above reduce to c-number factors. Indeed, r2 is a 
particular case of the general vertex function shown 
in Fig. [1] corresponding to n — 2, i. e. two of its legs 
are represented by single external lines (one for the 
constituent fermion and one for the constituent bo- 
son), while the third leg, for the physical fermion, is 
shown by a double line. Radiative corrections to each 
of the two external single lines are given by insertions 
of self-energy parts with their subsequent summation. 
The latter, of course, can be done directly within LFD 
by using the graph techniques rules, but we will choose 
a simpler way. 

Each on-energy-shell amplitude calculated in LFD 
must coincide with that found in the standard four- 
dimensional Feynman approach and taken on the mass 
shell. Hence, r2(s2 = "m?) coincides with its on-mass- 
shell Feynman counterpart ^. 

The summation of radiative corrections to external 
lines in the Feynman approach is technically easier 
than in LFD, since it can be done for each of the two 
lines independently. We allow the external particle 
momenta being off the mass shell (in order to avoid 
intermediate singularities), then sum up the radiative 
corrections and finally perform a limiting transition to 
the mass shell. Thus, summing a chain series of self- 
energy blocks - together with the mass-counterterm 
insertion - on the constituent fermionic line with the 
four-momentum ki brings the factor 



^ A three-leg vertex which enters, as an internal sub-block, in 
physically observed amplitudes is always off-shell. Taking it 
on shell, we imply its analytical continuation into a nonphys- 
ical kinematical region 
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X 



^1 



= lim 



— m 



1 - 



(17) 



where I]r(/^i) = ^(/i^i) ^ S(m). This factor is nothing 
else than Zf given by Eq. ( [T6| ). Analogous procedure 
for the boson line leads to the factor Zf,. The total 
factor which appears due to the dressing of the two 
constituent lines is therefore ZfZ},. 

Concerning the double fermionic line in its 
renormalization factor is related to the normalization 
condition for the state vector. The eigenvalue equa- 



tion ( 13 1 transforms into a homogeneous system of lin- 
ear integral equations for the vertex functions. Hence, 
the solution is determined up to an arbitrary common 
factor. In practice, for solving this system of equa- 
tions, it is convenient to fix the (constant) one-body 
Fock component 0i , so that the other components (r2 , 
Fa, etc.) become proportional to it. Then the double 
line in F2 brings the factor (pi which, in its turn, is 
determined by the normalization condition ([5| for the 
state vector. Since 0^ is just the norm of the one-body 
Fock sector, is equal to where /i is the first 
term in the sum (|6|. 

The relation between F2 and r2 becomes therefore 
on the energy shell 



(18) 



F2(s2 = m ) = V^i r2(s2 = m ) ZfZb 



Excluding f2(s2 = m^) from Eqs. (151 and (18), we 
find that the renormalization condition reads 



F2(S2 



(19) 



where F2 is expressed from the eigenvalue equa- 
tion (13 1, through the BCC 50- A similar discussion 



of the renormalization condition in terms of the one- 
body component was already done in Ref . [TD] . 

When we neglect the boson dressing by fermion- 
antifermion fluctuations, as we do in this work, the 



condition ( 19 ) finally reduces to 



F2(s2 =m ) ^ g^Ii 



(20) 



Note that there is a simple relation between the one- 
body normalization factor Ii and the field strength 
renormalization factor Zf: 



as shown in Appendix [A] 



C. Renormalization scheme 

The above conditions imposed on the BCC and 
MC are necessary in order to express physical observ- 
ables, like the electromagnetic form factors, through 
the measurable coupling constant and masses. As a 
consequence, one should expect full cancellation of di- 
vergences. 

Such a program could be realized in perturbation 
theory or nonperturbatively if the Fock space is not 
truncated. The latter case is hardly achieved in prac- 
tice. Usually, Fock space is truncated to a finite or- 
der N of admitted Fock sectors, and the cancellation 
of divergences is not anymore guaranteed. For in- 
stance, looking at Fig. [5] for the calculation of the 
fermion propagator in the second order of perturba- 
tion theory, one immediately realizes that the cancel- 
lation of divergences between the self-energy contribu- 
tion (of order two in the Fock decomposition) and the 
fermion MC (of order one) involves two different Fock 
sectors [2 . This means that, as a necessary condi- 



Sm 



FIG. 2: Renormalization of the fermion propagator in the 
second order of perturbation theory. 



tion for the cancellation of divergences, any MC and, 
more generally, any BCC should be associated with 
the number of particles present (or "in flight") in a 
given Fock sector. In other words, all MC's and BCC's 
must depend on the Fock sector under consideration. 
The original MC, 6m, and the fermion-boson BCC, 
go, should thus be extended each to a whole series: 



(21) 



90 
6m 



901 , 
6mi, 



(22a) 
(22b) 
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with I = l,2,...N. The quantities goi and 5mi are 
calculated by solving the systems of equations for the 
vertex functions in the N = 1, N ~ 2, N = 3, 
... approximations successively. This FSDR scheme 
has been proposed initially in Ref. [IIJ and devel- 
oped as a full renormalization scheme in Ref. An 
alternative approach, also in the Pauli-Villars (PV) 
regulated Yukawa model with the two-boson trunca- 
tion, but with a sector-independent renormalization 
scheme, was developed in Ref. |12) . 



Note that the series (22) does not imply that we 



have an infinite number of counterterms or bare pa- 
rameters. We still have the original ones go and Sm 
in the Hamiltonian we start with, but they have dif- 
ferent values according to the level of approximation 
used in the calculation. In the limit of an infinite N, 
and if the Fock sector expansion converges, goN and 
Smpf should turn to the true BCC and the MC, re- 
spectively. This is completely analogous to the case 
of perturbation theory where, at each order n, one 

determines gg"'' and Sm^"\ 

Apart from the mass and vertex radiative correc- 
tions, the third type of divergences arises from the 
field renormalization, i. e. from the constants Z f and 
Zf). The values of these constants should also depend 
on the maximal number N of particles kept in a given 
truncation. Consider for instance the vertex function 
r2 represented by Fig. [T] for n — 2. The dressing of 
the physical fermion leg (the factor yH^) should be 
calculated for the truncation to the A^-th order. The 
situation changes however for the constituent (sin- 
gle) fermion line. The state in which the constituent 
fermion is considered already contains one constituent 
boson. Hence, even if the boson line is not dressed, the 
dressing of the constituent fermion leg involves radia- 
tive corrections of order {N — 1). In other words, the 
dressing factor Zf for the constituent fermion leg must 
be calculated for the lower, (A'' — l)-body truncation. 
Otherwise, we would go beyond our approximation, 
since the effective number of particles in which the 
physical fermion can fluctuate would exceed N. 



Taking this into account, the relations ( 15 ) and ( 18 ) 



for a finite order truncation N (and in the absence of 
boson dressing) obtain the following form: 



4^^f2(,S2 = m2)^zf-i)=5, (23a) 



T,{s2 = m')^^l[''%is2^m')zf-'\ 

(23b) 

The superscripts (A^) and (A^ — 1) here and below 



just indicate the order of the Fock space truncation in 
which the corresponding quantities are calculated. 
It follows from Eqs. (23a) and (21) that the renor- 



malization condition ( 20 1 simply writes 



(24) 



in the absence of boson dressing. 

For the simplest case of the two-body truncation, 
N — 2, one thus gets 



Ty{s2^m } ^ g, 



(25) 



since — 1. We recover here the condition given 
in Ref. This condition is however valid only for 
iV = 2. 



III. YUKAWA MODEL IN THREE-BODY 
TRUNCATED FOCK SPACE 

A. Eigenvalue equations 

We consider in this study the Yukawa model: a spin- 
1/2 fermion interacting with massive spinless bosons. 
The regularization is provided by the PV method. In 
addition to physical particles, we introduce therefore 
one PV fermion and one PV boson with (large) masses 
mi and fii, respectively. This amounts to extend the 
physical Fock space to embrace negatively normal- 
ized PV particles [5]. The interaction Hamiltonian 
in Eq. ( 13 1 is given by 



with 



^jj' — -\- Tp 



PV, 



(26) 



where ip and ip are the free physical fermion and boson 
field operators, while ^ppv and ifpv are their PV part- 
ners, with negative norm. The bosons are supposed 
to be neutral. 

In the three-body truncation, the system of cou- 
pled equations for the vertex functions, derived from 
the eigenvalue equation ( |13[ ), is shown graphically in 
Fig. [3j On the r.-h. s. of the last equation, the sum 
of the diagrams with permutated boson legs appears, 
reflecting the symmetrization of the amplitude due to 
the identity of bosons. Expressing Fa through r2 by 
means of this equation, and substituting the result 
into the second equation, we can exclude the highest 
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FIG. 3: System of equations for the vertex functions in 
the case of three-body Fock space truncation. 




Analytically, these equations read 
u{pu)T\u{p) 



u{pii) {Vi 
u{ku) (1/3 



V2)u{p), (27a) 
Vi^ + V^)u{p), 

(27b) 



where the indices i and j refer to whether the line of a 
constituent fermion [i) or a constituent boson (j) cor- 
responds to a physical particle = 0) or to a PV 
one (i, j = 1). The term V45 means the sum of V4 and 
V5. The explicit expressions for Vi_6 are given in Ap- 
pendix[C| Note that the first equation in Fig. [3]is just 
a constraint which determines 5m^. The contribution 
V5 on the r.-h. s. of Fig. |4] involves the two-body self- 
energy depending, due to the departure off the energy 
shell, on a four-momentum k with fc^ ^ m^.This self- 
energy is decomposed, in CLFD, according to [5]: 



S(/^)=5o2 



A{k'^)+B{k^ 



C{k 



u)-k 



(28) 



where the factors m are here for convenience only. The 
coefficients B, and C are calculated in Appendix [B] 
In contrast to the two-body case, the system of 
equations for the three-body truncation is rather non- 
trivial. For example, its iteration generates all the 
graphs for the self-energy which contain one fermion 
and two bosons, including overlapping self-energy 
type diagrams. The number of such irreducible graphs 
is infinite. Some of them are shown in Fig. [5] The so- 
lution of the system of coupled equations incorporates 
the sum of these contributions to all orders. 
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FIG. 4: Reduced equation for the two-body vertex func- 
tion, obtained from that shown in Fig.[3]after the exclusion 
of the three-body component. 



FIG. 5: Radiative corrections to the self-energy. 



order vertex function Fa from the full system of equa- 
tions. We thus obtain a reduced equation for the two- 
body vertex function. Together with the first equation 
in Fig. |3j it forms a system of equations involving the 
vertex functions Fi and F2 only, as shown in Fig. |4] 



Note that due to the covariance of our approach, we 
can identify the contribution ^ (/j to the self-energy 
which explicitly depends on the light front plane ori- 
entation. If not regularized, the coefficient C{k'^) is 
quadratically divergent and needs a priori both PV 
fermion and boson regularization. After this regular- 
ization however, C(fc^) = for any values of the PV 
fermion and boson masses. This makes the two-body 
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self-energy identical to the result obtained in pertur- 
bation theory in the Feynman approach [7]. The con- 
tributions A and B do respect chiral symmetry in the 
sense that they are equal to zero when the constituent 
mass m, as well as the physical mass M (which, in our 
case, coincides with m), goes to zero, without the need 
of an extra PV boson. This is at variance with the 
standard formulation of LFD where it is claimed that 
an additional PV boson is needed, if the PV fermion 
mass is kept finite [T5] . 

The parameters 1702 and 5rn2 are taken from the 
N ^ 2 calculation [5]. They are given by 



2 

902 



Sm2 



9 



1 - ffV2 



[A{r 



where 



(29a) 
(29b) 

(30) 



with zq = 2to [A'im^) + B'im^)] . The norms of the 
one- and two-body Fock sectors, entering the normal- 
ization condition ([6]), are 



7^(2) _ 1 „2 T 



I. 



:2) 



(31a) 
(31b) 



In the two-body approximation mentioned above 
and discussed in detail in Ref. 2 , the two-body ver- 
tex function is automatically independent of uj since 
C(fc^) = 0. Moreover, it is a constant (i. e. it does not 
depend on the momenta of the constituent). Due to 



all these, the renormalization condition (25 1 directly 



leads to the relation ( 29a I between the bare and phys- 



ical coupling constants. In principle, nothing prevents 
r2 to be w-dependent, since it is an off-shell object, 
but this dependence must completely disappear on 
the energy shell, i.e. for S2 = rn? ■ It would be in- 
deed so if no Fock space truncation occurs. The lat- 
ter, in approximations higher than the two-body one 
(i. e. for TV = 3, 4, . . .), may cause some cj-dependence 
of r2 even on the energy shell, which immediately 
makes the general renormalization condition ( 19 1 am- 



biguous. If so, one has to insert into the light-front 
interaction Hamiltonian new counterterms which ex- 
plicitly depend on uj and cancels the w-dependence 
of r2(s2 = rn?)- Its explicit form will be given in 
the next subsection. Note that the explicit covari- 
ance of CLFD allows to separate the terms which de- 
pend on the light front plane orientation (i.e. on uj) 



from other contributions and establish the structure 
of these counterterms. This is not possible in ordinary 
LFD. 



B. Calculation of the two-body Fock component 

The method of solution is similar to that used in 
the calculation 2^ ioY N — 2. We first decompose the 
vertex functions in invariant amplitudes. The vertex 
functions on the l.-h. s.'s of Eqs. (27), being matrices 



in the spin indices, can be decomposed in a full set of 
spin matrices. This decomposition is very simple in 
CLFD and takes the form 



{mi~M^)^\u{pu)u{p), 



u{kii) 



u-p 



y-ip), 



(32) 



where Tp\ is a constant, and are invariant functions 
of particle momenta, with toq = and fiQ = /j. We 
denote temporarily, in the first of the above equations, 
the physical fermion mass by M in order to avoid sin- 
gularities, since the equations contain the combination 
r\/ (mf — M^) which becomes indeterminate for i = 
at M = m. Using M ^ m allows to take a smooth 
limit limM^m[ri/(mf — M^)] — ipl; after that one 
may set M ^ m. 

Each of the functions 65^2 depends on two invariant 
kinematical variables. As usual, we define a pair of 
variables, consisting of the longitudinal and transverse 
(with respect to the three-vector a; ) momenta: 



UJ-k2 

oj-p 



, Rj. = k2j. - xp±, 



(33) 



where ^2 is the boson four-momentum. Then are 
functions of x and i?^. 

The renormalization condition (24), for = 3, im- 
plies two conditions 



rn?) 




(34a) 
(34b) 



for the spin components of r2 at S2 ~ rn^ , where 
the two-body invariant energy squared S2 is expressed 
through and x as follows: 



S2 



Rl 



Rl 



(35) 
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One should emphasize that the renormaUzation con- 
ditions are imposed on the two-body vertex function 
corresponding to both physical constituents. The 



condition (34a) defines unambiguously gos- The con- 
dition ( 34b I is not verified automatically if the Fock 



space is truncated for TV > 3, unlike the case N = 2. 
We should thus enforce it by introducing an appropri- 
ate counterterm, as explained above (see also Ref . [2 ) . 
It corresponds to the following additional structure in 
the interaction Hamiltonian: 



6nr. 



(36) 



where is a constant adjusted to make Eq. (34b) 
true. The operator (/j/{iuj-d), in momentum space, 
leads to the appearance of a new three-leg vertex 



^/{u!-k) on each fermion-boson vertex with total in- 
coming momentum k. In principle, a similar new uj- 
dependent counterterm should be also added to the 
Hamiltonian in order to cancel the oj-dependence of 
Sm^, in full analogy with the cancellation of 6™(s2 — 
m?) [7 . However, as we will see below, 5m^ is needed 
only for a calculation in the four-body Fock space 
truncation. For this reason, working within the three- 
body truncation only, we may not bother about addi- 
tional counterterms excepting that given by Eq. (36). 



To solve the system of equations (27), we substi- 



tute the decompositions (32) into the expressions for 
Vi-Q given in Appendix |C[ then multiply Eqs. (27a) 
and (27b) by u{pii) and u(fcii), respectively, to the 



left and each of them by u{p) to the right, and sum 
over spin projections. We thus get 



J 



UJ ■ p 



{A^+m^){Vl+V2){r^+m), 

iA^ + m^)iV3 + V45 + Ve){^ + m). 



(37a) 
(37b) 



The system of matrix equations ( 37a ) and ( 37b ) can means of the relation 



be transformed into a homogeneous system of ten lin- 
ear integral equations for ten unknown functions (two 
ipl , four b''^ , and four bt^ ) ■ These equations are ob- 



tained by taking the trace of Eqs. (37a) and 



37b 



(six equations) , and by taking the trace of Eq. ( 37b 
after the multiplication of its both sides by (four 
equations) . 

In order to achieve the limit toi — > oo, it is conve- 
nient to replace the functions tpl and hi 2 by the new 

functions a^, /i^, and Hj according to 



mi{m + mi) ' 
rrii Hi — {l — X 



bY = —H, 

m 



2(1 



(38) 



A careful analysis shows that in this limit the PV mass 
TOi disappears from the equations written in terms of 
Qij, /i^, and Hf. These functions have therefore a finite 
limit. Below we will imply that the limit mi 00 is 
taken and a^, h^, and Hi denote the limiting values. 

For further simplification of the equations, it is con- 
venient to introduce new functions /ig 1 ^.nd H^ ^ by 



"0.1 

H^ 
-"0,1 




with 



1 



« = .903 - 



Using the substitution ( 39 ) and denoting 



Z' 



2Z^ 

.903 



(39) 



(40) 



(41) 



the initial system of ten equations splits into two sub- 
systems. The first one contains two equations involv- 
ing the ratio ai/ao. The value of Sm^ just ensures 
that both equations define the same quantity ai/ao. 
It is not interesting for our study in the three-body 
approximation, since ai, as will be seen below, drops 
out from the observables we calculate here, while ao 
is uniquely determined by the normalization condition 
for the state vector. As already mentioned, we also do 
not need to calculate Sm^ itself. It is used as an input 
in the calculation at the next, A'^ = 4, truncation. 
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The second subsystem of eight equations involves 
the eight functions hi and Hi^ only, since the ratio 
ai/ao is absorbed into the definition of Z'^ in Eq. (41|. 
We thus get 



K = 1 



/ [Kihi + Kfh{)+g''^l 



g'^ (-Kihi + Kih{)+g'\i, 



Hi = Zi{l-x) + 2-x 
+ 9'" {Ki 



(42) 



Hi 



1 



-KiHl+KlHi 



where 



and 



Ki 



KI = 



Ki - 



(1 + 5^^0)^ 



1 r 2[A(si)+S,(si)]m2 

— <^ Br[Si) ^ 



9 ' 

— si 

gr(si) 

m 



The substracted self-energy contributions Ar{si) and 
Sr(si), are given by 

Ar{si) = yi(si) - A{m^), Br{si) = 6(si) - B{m^) 



with 



Sl 



+ (1 - x)rr? 



(44) 



The functions ,4 and B are given in Appendix [Bj 
while the integral terms Jq and /g ^ are given in Ap- 
pendix [Dj 

The limit of infinite PV mass /ii is not easy to 
perform analytically, as it was done for mi. Set- 



ting /ii — ^ cx) directly in Eqs. (42) makes some in- 



tegration kernels singular (they decrease too slowly at 
^ oo). The dependence of physical observables, 
like the AMM, on /ii will therefore be studied numer- 
ically. 



Note that although g'^2 in Eq. ( 29a ) can become infi- 
nite (for J2 = 1/5^) and changes sign from positive to 
negative at sufficiently large values of the PV boson 
mass /ii , the eigenvalue equations ( 42 ) do not show 



any singularity when iJqj goes to infinity. Indeed, 

These equa- 
liil, with zo 



does not appear in the equations (42 1, 
tions depend only on g'^ , given by Eq. 
strictly positive. Therefore, g'^ is strictly positive and 
finite. We shall come back in Sec. |V]to the interpre- 
tation of the limit of large /ii, when both g^2 ^-nd the 

norm of the one-body sector I\ are negative, while 

(2) 

the norm of the two-body sector is larger than 1, 
from Eqs. 

The constant Z'^ entering the system of equa- 



tions ( 42 1 is determined from the renormalization con- 



dition (34b), while gga needed to calculate the renor- 
(43) malized vertex functions in Eqs. (39) is determined 



from the renormalization condition (34a 



The components bf^2i.^2 — w^) entering these renor- 
malization conditions are expressed through the so- 
lution of the system of equations (42) by means of 
Eqs. (38 1 and (39). The kinematical point S2 — rri^ 
belongs to a nonphysical region, but there is no need 
to make an analytical continuation to this region of 
the solution and Hi^ found numerically. Indeed, 
the integral terms in Eqs. (42) involve integrations 



within the physical domain only. One can simply set 
j — 0, R± = R*_i, X = X* , where R*^ and x* are 
determined by the condition S2 = rn^, and calculate 
the integral terms by substituting there the previously 
found solution hl{R±,x) and Hj{R±,x) for physical 
values Rj_ and x. After that Eqs. ( [42| reduces to 
a system of four ordinary linear inhomogeneous equa- 
tions for h^{R*j_,x*) and H°iR*j_,x*). Finally, relating 
the calculated quantities hQ{R*j_,x*) a nd H n(R* , x*) 
to 6?°2(*2 — ™^) we get goa from Eq. (34a)'^ and Z'^ 
from Eq. (Isibl). 



The condition S2 — rn however does not determine 
and X* simultaneously. It is convenient to fix x* 
from Eq. (|35|: 



R*j 

somehow and then find R 

R*l = - [x*'^m 



(1 



X )fl 



(45) 



Since the two-body vertex function ( 32 ) on the energy 
shell must turn into a constant, the functions 6i'2(s2 — 



^ More precisely, we get not go3 alone but the product gosctQ. 
The quantity ao is found from the normalization condition 
for the state vector. This procedure requires knowing the 
three-body normalization integral which is calculated in the 
next section. 
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m^) also must be constants. In other words, if one spin structures as follows: 



relates the arguments of these functions by Eq. ( 45 1 , 
their values are independent of the choice of x* . It 
would be so in exact calculations, i. e. if Fock space 
was not truncated. A finite order truncation makes 
the Fock components, even at S2 — rn^, a;*-dependent. 
As advocated in Ref. [2j, we choose x* = We 

shall investigate in Sec. [v|how 6™2(s2 = m'^) depends 
on the choice of x* . 



C. Representation of the three-body component 

We can find the three-body component by calculat- 
ing the amplitude corresponding to the r.-h. s. of the 
equation shown by the last line in Fig. [3j 

The general form of the relativistic vertex function 
of a system composed from one constituent fermion 
and two spinless bosons with total spin 1/2 reads 



<(A:i)r„^(l,2,3)«^(p), 



(46) 



where rQ,^(l, 2,3) is a 4 x 4-matrix in the indices a, /3. 
The arguments of Fa, denoted symbolically by num- 
bers, mean three pairs of the standard variables 



Rz_L = kj_L - XiPj_, 



uj-p 



with I = 1 corresponding to the fermion and / = 2, 3 
to bosons. Here u'^_^{ki) is the bispinor of the con- 
stituent fermion, (p) is the bispinor of the physical 
fermion (of the composite system), ai, a are their spin 
projections in the corresponding rest frame. Since 
CTi = ±1/2 and u — ±1/2, we have in general 2x2 = 4 
matrix elements. Usually, parity conservation reduces 
this number by a factor of two. However, this is not 
the case in relativistic calculations, for a n-body wave 
function with n > 3 [14j . This wave function is deter- 
mined by jour independent matrix elements or, equiv- 
alently, by four scalar functions like 

u(fci)F3(l,2,3)w(p) (47) 
= u{ki){g\ Sx +52 'S'2 +53 S-i ^ g^S^u{p). 

For simplicity, we omitted the bispinor indices and the 
indices marking the particle type (cither physical or 
PV one). It is convenient to construct the four basis 





= 2x1 - 




+ xim 


ijjp 








S2 


= m — 

UJp 


1 






S3 




2x1 


- {rui ~ 


xim) 


S4 


— imC 




75 





uj-p 



75, 



UJp 



(48) 



with Xl — ^^^^ and m,- being the internal fermion 
mass (either physical or PV one, depending on which 
type of fermion the momentum ki corresponds to), 
while Cps is the following pseudoscalar: 



Cps — 



1 



rn^ujp 



''hflhl^Ppi^'y- 



(49) 



The function Cps can only be constructed with four 
independent four- vectors. This is the case in LFD for 
n > 3. In the non-relativistic limit, one would need 
n > 4. We can then construct two additional spin 
structures S3 and S4 of the same parity as Si and 
S'2 by combining Cps with parity negative matrices 
constructed from Si, S'2, and 75 matrices. Instead 
of k2p,k3iy one could have taken any pair of momenta 
[kifik^i, or ki^k2v). We take the boson momenta for 
symmetry. With this definition 



1 



ps 



[^2_L^3_L - (R2_L-R3_l)^]- 



(50) 



The three-body vertex function r3(l,2, 3) is com- 
pletely determined by the four scalar functions 
31-4(1,2,3) in Eq. (47). They depend on Ri_3_j^ in 
the form of their scalar products among themselves 
and on xi_3. Since 

RlJ. + R2_L + R3_L = 0, Xi+X2 + X3 = l, (51) 

we can exclude, for instance, Ri_l and xi. The func- 
tions gi^2 are symmetric relative to the permutation 
203, whereas 33^4 are antisymmetric: 

51,2(1,2,3) = 51.2(1,3,2), 
53,4(1,2,3) 53.4(1,3,2), 

so that the product Cps 53,4(1, 2, 3) which appears in 
S3 53 and S4 54 is symmetric. 

Each component 5„ is represented as a sum or a 
difference of two terms: 

51,2(1,2,3) =51,2(1, 2, 3) +51,2(1, 3, 2), 
53,4(1, 2, 3) = 53,4(1, 2, 3) - 53,4(1, 3, 2), (52) 
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where the permutation 2^3 means 

In their turn, g„(l,2, 3), according to the last hne in 
Fig. [Sj are hnearly expressed through the functions 
/iq^j^ and Hq\ which form a sohition of the equations 

5«(1,2,3) = ao K go2 a„o(l, 2, 3)/i^^ (2) 

+a„i(l, 2, 3)/if (2) + A„o(l, 2, 3)i?^'^ (2) 

+ A„i(l,2,3)#f (2)] . (53) 

The coefficients a and A in this formula are given in 
Appendix |Ej 

We can finally calculate the three-body normaliza- 
tion integral. It is given by 



h = 



-.3233 



32:3-3=0 



(S3 — W?)' 



with 



-Tr[f3(^i+m)r3(:^ + m)] 



dD:, (54) 



(55) 



and, as usual, F = 7°F^^7°. The factor ^ in Eq. (55) re- 



sults from averaging over ini tial state spin projections, 
while the factor ^ in Eq. (54 1 is the combinatorial 
factor (^n-iy. oi'igiiiating from the identity of the two 
bosons. The contribution of PV fermion is omitted 



since it disappears in the limit mi — > oo. The phase 
space volume element has the form (see Eq. (3.23) 
from Ref. 0): 



R-2± 

3 



R 



(56) 



X 5{xi + X2 ^ x^i - 



(PRi±dxi 



IV. ELECTROMAGNETIC FORM FACTORS 

A. Electromagnetic vertex in CLFD 

The electromagnetic vertex contains contributions 
of one-, two-, and three-body Fock sectors, as shown 
in Fig. [6j They are expressed, in our FSDR scheme, 
in terms of the external electromagnetic BCC cqi, as 
explained in Ref. [2] . These coupling constants are all 
identical to the physical fermion charge, i. e., Cqi = e 
for all Vs. Note that this important property of QED 
is not preserved in general if FSDR is not used. 




FIG. 6: One-, two-, and three-body contributions to the 
electromagnetic vertex. 

The decomposition of the spin- 1/2 electromagnetic 
vertex in CLFD is given by [SI [E] 



u{p')GPu{p) = eu{p') 



2m \^'P 



uj-p 



B: 



u{p) (57) 



with P — p + p' , and q — p' — p. Fi and F2 are the 
physical form factors, while i3i,2,3 are spurious (non- 
physical) contributions which appear if rotational in- 
variance is broken, e. g. by Fock space truncation. 



The decomposition (57) enables to separate unam- 



biguously the physical form factors from the nonphys- 
ical ones. Under the condition oj-q = 0, all F1.2, S1-3 



depend on Q'^ 



-q only. 



We shall represent q = {qo,A,q\\), where and 



A are, respectively, the longitudinal and transverse 
components of the momentum transfer with respect 
to the three- vector uj. Since u-q = 1^0(90 ~ ^n) = Oi 
have = A^. 

After construction of the matrix 



4to^ 



m)GP{^ + m), 



(58) 
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and calculation of the traces 
C4 = Ti[OpUj'^]m/uj-p, C5 = Tr[OpCj''(Z)]m^/(cj-p)^, 



the electromagnetic form factors write [7] 



(59) 



1 2m^ 

Fi = -C5, F2 = ^(C5-C4). (60) 

The value Fi{Q^ = 0) equals one, since it coincides 
with the norm of the state vector. The value F2{Q^ — 
0) is just the AMM. 



B. One-body contribution 



form factors in terms of the solutions /i^, fff of the 



system of eigenvalue equations ( 42 ) . After taking the 



limit mi — )■ 00 the result is as follows: 



2b 



9 o 1 



{-ly / d'^Ri_ / dx 



[R\ + x2m2 + (1 - + x^m^ + (1 - x)^i'^] 



(63) 



^2.26 — 



47r3A2 



^(-1)^ / d^Ri_ f dx 



The one-body contribution to the form factor Fi is 
given by the first diagram in Fig.|6] It does not depend 
on Q2 and coincides with the norm of the one-body 
sector: 



(61) 



There is no one-body contribution to the form factor 

F2- 



F2.U = 0. 



(62) 



[Rl + x2m2 + (1 - x)n^^][R'i + x^m^ + (1 - x)^2] 

(64) 

Functions with primes depend on R' ^ and x. The 
value Fi^2b{Q^ — 0) coincides with the two-body con- 
tribution to the normalization integral, 

^2^2 1 fOO fl 

r.~^^{-iy / d'R^ / dxx 



167r3 



C. Two-body contribution 

The two-body contribution to the electromagnetic 
vertex, as given by the second diagram in Fig.|6j writes 



1 

(27r)3 
d^Rj_ 



E (-1: 



i+i'+j 



i.i' .j—0 
r-1 



dx 



/o 2a;(l-x)2 



/i'j 
2 



where T2 is given by Eq. |32| with b'l^ = ^1^2 



kii {k[^i) is the momentum of the constituent fermion 
incoming to (outgoing from) the elementary electro- 
magnetic vertex, S2 = {kii+k2j)^ , S2 = {k[^,+k2j)'^ , 
and k2j is the constituent boson momentum. T2 has 
the same decomposition as T2, with the replacement 
61^2(^1,2;) b(^2iR'l^x) with = R_L - xA. 



Using the relations (38 1 and (391, we can calcu- 



late the two-body contribution to the electromagnetic 



R 



I (hi) 



-"0 



[R\ + a;2TO2 + (1 



a:)M|]2- 



(65) 



To calculate the two-body contribution to the 
AMM, which is given by the value J2,2fc(<3^ — 0), one 
should go over to the limit A — > in Eq. ( 64 ). The cor- 



responding analytic formula includes derivatives over 
R±^ from the Fock components. For numerical calcu- 
lations it is however more convenient to find -^2,26 at 
small but finite and then, decreasing the latter, to 
reach desired accuracy. The result of this numerical 
limiting procedure is very stable. 



D. Three-body contribution 

The three-body contribution to the electromagnetic 
vertex reads 

u{p')Gl,u{p) (-1)'^^'^ 

u{p')f'zm[ + m)Y{h + m)V^u{p) 



4x2x2X3 (s3 — m?){s'^ — m?) 



dD. 
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where dD^ is defined by Eq. ( 56 1 and 

m2 Rl^+^^l ^3±+m: 



S3 



rIj_ 



X2 



X3 



(66) 



S3 differs from S3 by tlie following shift of the argu- 
ments: 



R-i_L - 


^ R' 


i± 


= R.i± 


+ (l-a;i)A, 




R-2_L - 


^ R' 


2± 




-a;2A, 




R-3_L - 


^ R' 


3± 


= R-3± 


-a;3A. 


(67) 



r'3 is obtained from r3 by the same shift of argu- 
ments. From the decomposition (|47|, we can calculate 



C, and construct the matrix by means of Eq. ( 58 1 



The form factors are thus given by Eqs. (60) and read 



Fi,3b = / iC[\^9i9[ + Cg^9292 + 5333 



Cil^9m + 2Cg^939'i) 



di 



(68a) 



2,36 



2 1 {C[f9i92 + Cif9A9'i + Cg^9392 



C349394 



dDs 



(68b) 



where 



di = ■m'^X2Xz{m'^ — S3){m^ — S3), ^2 ~ 



2A2 



gl^ differs from gn by the shift of the argum ents ( 67 1 . 

(1 2) J L J I 

The coefficients in Eqs. (68a I and (68b) are 

given in Appendix [F] 

The value Fi^-^biQ^ = 0) coincides with the norm, 
/3, of the three-body sector given by Eq. (54 1. The 



quantity ao which has been unknown up to now, is 
determined from the normalization condition for the 
state vector: 



NUMERICAL RESULTS 



The solution of Eqs. (42 1 is found by a matrix inver- 



sion after discretization of the integrals, using Gaus- 
sian method. All integrals are finite at finite PV bo- 
son mass /ii. As already mentioned, the limit of infi- 
nite PV fermion mass mi has been done analytically, 
while the Fock components /i^. Hi and, hence, 1)^2 in 
Eq. (32 1 do depend on the PV boson mass ^1. The 



numerical calculations have been performed on an or- 
dinary modern laptop. 

The AMM is calculated for a typical set of physical 

parameters m = 0.938 GeV, ^ = 0.138 GeV, and two 

2 

values of the coupling constant a = = 0.2 and 
0.5. This mimics, to some extent, a physical nucleon 
coupled to scalar "pions". The typical pion- nucleon 
coupling constant is given by i? = jp-{k) where (k) 
is a typical momentum scale, and gA and F^^ are the 
axial coupling constant and the pion decay constant, 
respectively. For (fc) = 0.2 GeV we just get a ~ 0.2. 



o 

uT 0,02 



log 



0,12 
0,10 

0,08 
^ 0,06 

o 

uT 0,04 
0,02 
0,00 



log [i^iVn^] 



"0 + ^2 + ^3 - 1- 



(69) 



Since both I2 and are proportional to Oq, then, de 
noting J2,3 = q;oK^/2,3, where k is defined by Eq. (40 1 
we immediately get 



1 



l+K^{l2+h) 



(70) 



FIG. 7: The anomalous magnetic moment in the Yukawa 
model as a function of the PV mass /xi, for two different 
values of the coupling constant, a = 0.2 (upper plot) and 
0.5 (lower plot). The dashed and dotted lines are, respec- 
tively, the two- and three-body contributions, while the 
solid line is the total result. The AMM value calculated 
in the N — 2 approximation is shown by the thin line on 
the right axis. 
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We 



We plot in Fig. jrjthe AMM as a function of log 

for the two different values of a pointed out above, 
show also on each of these plots the value of the AMM 
calculated for the N — 2 truncation, which coincides 
with the AMM obtained in the second order of per- 
turbation theory. The results for a — 0.2 show rather 
good convergence as /ii — > oo. The contribution of 
the three-body Fock sector to the AMM is sizeable 
but small, indicating that the Fock decomposition ^ 
converges rapidly. This may show that once higher 
Fock components are small, we can achieve practi- 
cally converging calculation of the AMM. Note that 
this value of a is not particularly small: it is about 30 
times the electromagnetic coupling, and is about the 
size of the typical pion-nucleon coupling in a nucleus. 

When a increases, we see that the contribution of 
the three-body sector considerably increases. For a — 
0.5 the three-body contribution to the AMM starts to 
dominate at large values of fii. The dependence of 
the AMM on the PV boson mass /Lti becomes more 
appreciable, although it keeps rather small. 



the contributions of the one-, two-, and three-body 
Fock sectors to the norm of the state vector for the 
two values of the coupling constant, considered in this 
work. We see again that at a = 0.2 the three-body 
contribution to the norm is small, while it is not neg- 
ligible, and increases with fii, when a — 0.5. 

At very large values of /ii, and for large a, Ii be- 
comes negative. As already mentioned, we still get 
a well defined solution of Eqs. (42), and there is no 



discontinuity whatsoever in the value of the AMM. 
As shown in Fig. [Tj the convergence of the AMM as 
a function of the PV boson mass is expected in any 
case to settle much before we enter into this regime. 
According to renormalization theory, the mass of the 
PV boson should be much larger than any intrinsic 
momentum scale present in the calculation of phys- 
ical observables. With this limitation, physical ob- 
servables should be independent of any variation of 
the PV boson mass, within an accuracy which can 
be increased at will. This is what we found in our 
numerical calculation for small enough values of a. 
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FIG. 8: Individual contributions of the one- (dashed line), 
two- (dotted line), and three-body (solid line) Fock sectors 
to the norm of the state vector as a function of the PV 
boson mass fii, for a — 0.2 (upper plot) and a = 0.5 
(lower plot). 

In order to have a more physical insight into the 
relative importance of different Fock sectors in the de- 
composition ^ for the state vector, we plot in Fig. [s] 
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FIG. 9: The spin component b°° 
function ( 32 1 calculated at S2 = 
for a = 0.2 (dashed line) and a 
typical value of /ii — 100 GeV. 



of the two-body vertex 
, as a function of x, 
= 0.5 (solid line), for a 



In order to understand the possible origin of the 
residual de pen dence of the AMM on /ii, we plot in 
Figs. [9] and 10 the two-body spin components 5™ and 

as a function of x. As 
(s2 — vn?) and 



III 



calculated at S2 
we already mentioned in Sec. 

^2*^(^2 = w^) should be independent of x in an exact 
calculation. Moreover, should be zero. It is here 

— by 



fixed to zero at a given value of x ~ x* 
the adjustment of the constant Z'^ in the system of 
equations (42 1. We clearly see in these figures that 
is not a constant, although its dependence on x is 
always weak, while 63° is not identically zero, although 
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FIG. 10: The same as in Fig. [9] but for the spin component 



FIG. 12: The same as in Fig. |11[ but for the component 



its value is relatively smaller than that of fej" ioi a 
0.2, and starts to be not negligible for a = 0.5. 
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FIG. 13: The component defined by Eq. (391, as a 
function of x at R± = 0, for a typical value of — 100 
GeV and for a = 0.2 (dashed line) and a = 0.5 (solid line). 



FIG. 11: The component h'^ defined by Eq. ( [39| , as a 
function of R± at 2: = ^^''^ ^ typical value of ni — 100 

GeV and for a = 0.2 (dashed line) and a = 0.5 (solid line). 



We plot in Figs. 11 and 12 the two physical compo- 
nents /in and HS as a function of i? 1 , at x = . As 



expected from the system of eigenvalue equations ( 42 1 



the functions Hq and Hq tend to constants at large 
R±. Hence, the functions and Hq related to them 



by Eq. (391 tend to constants too. Note that the two- 
body wave function (j)2 related to r2 by Eq. ( 11 ) goes 



to zero at large momenta due to the rapidly decreasing 
kinematical factor (s2 — m^)~^. 

For completness, we plot in Figs. [13] and [M] the two 
physical components /ig and Hq as a function of x, at 
= 0. 
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FIG. 14: The same as in Fig. 1131 but for the component 
-no- 



VI. CONCLUSION 



We calculated the anomalous magnetic moment of 
a fermion in the Yukawa model, in the first non-trivial 
approximation, incorporating a constituent fermion 
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coupled to zero, one, and two scalar bosons, i.e. within 
the three-body Fock state truncation. We applied a 
general formalism based on the covariant formulation 
of light-front dynamics and an appropriate Fock sector 
dependent renormalization scheme which enables us 
to control uncancelled divergences when Fock space is 
truncated. We paid particular attention to the renor- 
malization conditions necessary to relate the bare cou- 
pling constant to the physical one. To do this, we need 
to identify all spurious contributions originating from 
the violation of rotational invariance, coming from the 
Fock space truncation. This is possible in our covari- 
ant formulation. 

The anomalous magnetic moment shows a very nice 
convergence as a function of the regularization scale 
(the Pauli-Villars boson mass ^i), for the coupling 
constant value a — 0.2 which mimics a nucleon cou- 
pled to a scalar "pion". For this value of a, the two- 
body component gives a dominant contribution to the 
anomalous magnetic moment. As a increases, we see 
the onset of higher Fock components. 

This shows up in the large contribution of the 
three-body component, and in the dependence of the 
anomalous magnetic moment as a function of the reg- 
ularization scale. We believe that this latter depen- 
dence should be largely, if not completely, removed by 
incorporating the relevant fermion-antifermion contri- 
butions to the three-body Fock components. We are 
presently investigating these contributions [IB] . 

As we have seen in our study, the calculation of 
nonperturbative properties of bound state systems de- 
mand to control all approximations in a quantitative 
way, in order to be able to make physical predictions 
order by order in the Fock space truncation. We think 
that the combination of the covariant formulation of 
light-front dynamics with an appropriate Fock sector 
dependent renormalization scheme is a quite promis- 
ing method to investigate these properties in a very 
elegant way, with a minimum of Fock components and 
computational time. 
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Appendix A: Relation between the Fock 
component normalization and the field strength 
renormalization factors 



We shall prove here the relation (21) in the general 
case, i. e. without Fock space truncation. We omit 
for a moment antifermion contributions generated by 
the process b — ?> //. They will be incorporated be- 
low. The fcrmion self-energy is given by a sum of ir- 
reducible graphs with all possible intermediate states: 



n=l 



(Al) 



Here Snd*) is the contribution from a graph with 
n intermediate states. For example, the first graph 
in Fig. [5] for the self-energy contains one intermedi- 
ate state, the second one contains three intermediate 
states while the third one contains 11 intermediate 
states. 

For the calculations we will use the 3D light-front 
graph techniques exposed in Ref. [3j . It can be repre- 
sented in different equivalent forms. We use the form 
in which an amplitude is represented as a product of 
energetic denominators (one denominator for each in- 
termediate state) multiplied by appropriate spin ma- 
trices. The amplitude (with a conventional minus sign 
for the self-energy) corresponding to a graph with n 
intermediate states has the following form: 



'9o 



-jdD n 

i—l 



,2 



P 



(A2) 



Since we do not truncate Fock space, we use the same 
BCC in all vertices. Si is the square of the invariant 
energy of i-th intermediate state. The product Y[i=i 
runs through all n intermediate states. The integra- 
tion in Eq. ( A2 ) is performed over all independent 
variables. 

The decomposition of S(^) is similar to that given 
by Eq. p): 



S](^) =^(p2)+6(p2)i+c(p2) !^+Ci(p2)^, (A3) 



ui-p 



where 



Two of us (V. A. K. and A. V. S.) are sincerely 
grateful for the warm hospitality of the Laboratoire 
de Physique Corpusculaire, Universite Blaise Pascal, 
in Clermont-Ferrand, where part of the present study 
was performed. 



Ec^. ( A3 1 is the most general form of the fermion self- 
energy m CLFD. The term with the function Ci{p'^) 
does not appear for the two-body {N — 2) Fock space 
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truncation, but it may appear for > 3. We give 
here the coefficient C which will be used below: 



1 

Am 



-Tr 



P 

uj-p 



(A4) 



To find ^ 



we first replace in Eq. ( A3 ) p 



by M^, 1^ by M (that is, replace by p" /M) and 
then calculate the derivative over M . ft is convenient 
to make this replacement by using the formula 

s(M) = ^Tr [(^ + hmmp^^M^ ■ 



We get 



d i) 



(A5) 



We substitute here S(^) from Eqs. (Al) and (A2) 



The contribution of derivative from the j-th fac- 
of the denominator, which results from 



-A/2 



tor 



Eq. (A2), reads 



d i> 



= -go+^2M I dD Tr 



AM 



n 



(Sj - M2)2 



n 



(A6) 



The factor 



9i 



I n 



M2) 



(A7) 



is a contribution of the graph with j — 1 intermediate 
states into the vertex function, and similarly for the 
second product. In contrast to Eq. (A6), where the 



integration dD runs over the phase volumes of all the 
intermediate states, the integration dD' in Eq. (A7) 



runs over the phase volumes of the intermediate states 
«i = 1, • ■ • , j - 1 only. 

Since all the four-momenta are on the corresponding 



mass shells fc^ — m'j, we have 



a=±l/2 

^-Tr[i^ + M)0] - ^ E ^^ip)Oux{p), 

A=±l/2 



for an arbitrary matrix O. The factor ^ in the last 
equation is introduced for averaging over initial spin 
projections. 

We therefore get 



d i) 



i,=M 



2 ^ 



dD 



ux{p)V;jU„{kj) [ux{p)TjU„{kj)\^ 
{s, - A'P) {s, NP) 



2^ 

A.cr 



dD, (AS) 



Here the integration dDj runs over the phase volume 
of the j-th. intermediate state not included in the in- 
tegral for Tj . The vertex function Tj may correspond 
to any fixed number of particles in the state j allowed 
by a given graph. We took into account that the fac- 
tor l/(sj — A/2) turns each F into (p, according to 



Eq. (11). Taking the sum over all the graphs and over 
all the intermediate states j, we recover in Eq. (AS) 



the contribution to the normalization integral In>2 
from all the A^-body states with N > 2 (each inter- 
mediate state in irreducibles graphs for S contains at 
least two particles). Since the rules of the graph tech- 



19 



niques used to calculate E imply that the one-body 
states are normalized to 1, this means that In>2 cor- 
responds to a state vector normalized by the condition 
/i = 1. If /i ^ 1, then Eq. (A6| determines the ratio 
—In>2/Ii- That is 



The contribution of the derivative of other factors 

reads 



in Eq. (A5), except for JJ^ 



-An ' 



d i) 



lN>2 



(A9) 



i>=M 



dD 



dM 



{ ^Tv[{i, + M) Y^Uih + ^^)V=M^ ] 



(AlO) 



Consider first the case when the products in Eq. (AlO ) 
contain only one factor with a fixed i. Then E cor- 
responds to the first graph in Fig. [5] We calculate 
the trace in Eq. (AlO), using the following explicit 



expression for the scalar product: 



where the variables Rij_ and constructed ac- 

cording to Eq. (33). Calculating then the derivative 
over M, we find 



1 



d i) 



-9l 



dD- 



dM l4M 



Tr[(^+M)(/^, + m,)]} 



-9l 



dD- 



(All) 



where d-Do 



d Ri±dxi 



(27r)3 2xi(l- 

space volume element. 

Let us calculate now the value of the coefficient C 



is the two-body phase in Eq. (A3). It is given by Eq. (A4). We still consider 



a particular case and keep one factor only with fixed 
i. Then C obtains the form 



dD- 



4m 



dD- 



2xIm^ (-^J-l + 



(A12) 



Comparing Eq. (All) with Eq. (A12), we find the 
relation 



d i> 



i)=M 



C. 



It turns out that Eq. (A13) is valid in the most 



(A13) 



general case. In the latter case, but still without an- 



tiferm ions, we get in the numerator in Eqs. (All ) and 
(A12) a product of the matrices n"=i(/^i + in- 
stead of the single term {^^ + rrii). This product can 
be decomposed in the full set of the 4x4 matrices as 
follows: 
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i=l 



(AM) 



where aih^.h^) = ^{fi.^fi^^ - The coeffi- 

cients Gi-a depend on the scalar products of the four- 
momenta fci, . . . fc„: 



1 



\Xi ra. 



+ (xjRij_ 



a;iRjj_)^ -I- x^.m}] 



It is important that these scalar products and, hence, 
the functions Gi_3 do not depend on M . Therefore 
can be extracted from the opera tor -^ jj ■ 

We 

replace {^^ + rrii) in Eqs. (All) and (Al2i b y the 
product ri"=i(/^i +'mi) represented in the form ( A14) 
(and take the product of the denominators). 'I'he 
matrices 75 and 75 ji^ give zero contributions to both 
Eqs. (All) and (A12), whereas with the matrices 1, 
/ifj, and a{ki^,ki^) we reproduce the relation ( |A13 1. 

The incorporation of antifermions (say, the /// 
intermediate state, in addition to bosons) does not 
change the form of the denominator (though the en- 
ergies Si incorporate now the antifermion momenta). 
That results in Eq. ( A9 ) . The corresponding numera- 



tor contains now the spin matrices of all the fermions 
-I- antifermions [we get one factor \/{si — p^)_and a 
product of three matrices (± l^^+m.i) for the /// state; 
the signs plus and minus stand for fermions and an- 
tifermions, respectively]. We still can decompose the 
full numerator according to Eq. (A14) and again re- 



produce the formula (A13) 



In this way, taking the sum of Eqs. (A9) and (A13| 
we finally find 



d i> 



LN>2 



h 



A/2 



C. 



(A15) 



If rotational invariance is preserved (it can be violated 
for instance by omitting some time-ordered graphs or 
by using rotationally non-invariant cutoffs), C is zero. 
It is indeed zero, for example, in the two-body ap- 



proximation with the PV regularization [see Eq. ( B6 ) 



in Appendix IB]. If C = 0, substituting Eq. (A15) into 



Eq. (16) and taking into account that Ii + In>2 ~ Ij 



we prove the relation (21) 



Appendix B: Self-energy coefficients 

We give here explicit formulas for the coefficients 
yl(fc2), S(fc2), and C(fc2) entering Eq. ^ for the 
two-body self-energy in the Yukawa model. If S(/{f) 
is known, these coefficients can be found as follows: 



9I2 C{k') = ^Tr 



Tr[S(^)(^] 



E(^) ^- i/;- 



i-k 



(Bl) 
(B2) 

.(B3) 



In the Yukawa model, the self-energy regularized by 
one PV boson and one PV fermion reads, in CLFD, 



S(« 



9I2 
(27r)3 



cPR\ 



dx 



Q 2a::(l — x) 

{j^l + TUi) 



(B4) 



where ki is the internal fermion four-momentum. The 
light-front variables are, as usual, Rj^ — \i2i_ — a;kj^, 
X — uj-k2/uj-k {k2 is the boson four-momentum), and 



1 



(B5) 



The coefficients A and B converge without PV 
fermion (i. e., they have finite limit when mi — > 00). 
Substituting Eq. ([bI]) into Eqs. (|B1]) and ([^, inte- 
grating over (PR± and omitting the PV fermion, we 
get 
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log 



167r2 



167r2 

1 

(1 - x) loi 





a;(l - x)fc2 + (1 - x)^!^ 
x{l - a;)fc2 + (1 - x)iil 
x{l - x)P + (1 - x)n'^ 
x{l - a;)fc2 + (1 - x)fil 



dx. 



dx. 



Notice that in the limit /ii — > oo and for fixed fc^ A similar calculation of C{k^) requires, for con- 

the values of Arik"^) = A{k^) — A{m'^) and Brik"^) = vergence, not only one PV boson, but also one PV 
B{k'^) — B{m^) are finite. fermion. We thus find 



32m7r2 /n l~xjn ' R"^. - x(l - x)k'^ + (I - x)n^, + xmf 



Appendix C: Right-hand sides of the eigenvalue of these equations are denoted by u(pii) (Vi + V2) 

equations ([27]) and M(fcn) (V3 + V45 + Vg) 'u(p), respectively. The ex- 

plicit form of Vi_6 is the following: 



The system of equations (27) determines the one- 
and two-body Fock components V\, V^^ . The r.-h. s.'s 



V, = ^rn3E(-ir' 2 (Cla) 

_ , v^. /" d'^R'i_ dx' {}([,, +mi>) ^,'j' 

V2 - .903 1) J J^J^ 2^ 2K^ ^ ' 

^3 - 5^3E(-i)^'^|^ri, (cib) 

4(0; • - - x~ x')Ti'j>Ti"jj' ^ 

(Cld) 

with (7q3 = f;o3 + Zrp ^^"-^ obvious notations for the The term V45 stands for the sum of the contributions 
momenta of the particles in the intermediate states. 
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of the graphs V4 and V5 in Fig. [5] The two-body vertex 
functions r2 inside the integrands depend on and 
x', while those which are not integrated depend on 
Rj_ and x. After calc ulati ng the traces taken from the 
equations (37a) and (37b), we obtain scalar products 



which are expressed through the variables R_l , x and 
R^, a;'. Examples are given in appendix C of Ref. [2]. 

The values of r's, which appear in the above formu- 
las, are related to the invariant energies in the corre- 
sponding intermediate states. For example, Ti'/jji in 



Eq. (Cld) for Vq has the form 



where 



Si"jj' 



(fci 



2uj-p 



«2j 



2j' 



kii" 



k2j, and k2j, are the four-momenta in the inter- 



mediate states while s, for any intermediate state, is 
expressed through the light-front variables as follows: 




where R^j^ and Xi are constructed according to 



Eq. ( 33 1 . They satisfy the conservation laws similar 



to Eq. (511 



Appendix D: The integral terms in the 
equations ( 42 ) 



The numerators and denominators of the kernels in 



the integrals in Eqs. (42) are linear functions of the 



scalar products Rj^-R'j^ = R±R'j^ cos (j)' , where R'j^ is 
the integration variable. We can therefore analytically 
integrate over d(j)' , using the formulas 



Jo - 
Ji = 



del)' _ sign(A) 

2nD{A + B cos </.') ~ DVA^ - ' 

2tt 



COS ( 



1 



A 



;o 2TrD{A + Bcosct)') DB B 
One should substitute here 



Jo- 



A = R']_{l-x)x + x'[x{x + x')m^ + R'i{l-x') 



- (x + x' - l){x'fl^j + X^i^j,), 



B = 2R'j^R±x'x, 
D = -87r2(l-a;'). 



With 



= m^x'^ + {I - x')^i], + R'l, 
the integral terms obtain the form 

ii = / / dx' {-ly' (Di) 

Jo ij,^o 

X (cmh{+CmH{^, 

poo nl — X 

li ^ / / dx' j2 i-^y' (D2) 

Jo Jo ij,^o 



for n — 0, 1. These integrals converge due to the PV 
regularization (the sum over j'). The sixteen coeffi- 
cients c, C, c', and C depend on j' . They are given 
below. 



Coo 



-Ir'^R^xx' Jq 
1 



Ji 



-R'\{x - l)x + x'[x{-x{x' - 3) 



coi 
Coo 

Cm 



2,x' ~ A)m^ + R\{1 -x')] 
{x — l){x' — l){x' + x^jL^j,) I 
R' 

--^x{2x + x' -2)Ji , 
xx'[JaR^^{?>x' — 2) 



Ri_rii 
i?l(2-3a:)Ji], 



[R'Ax-l)Ji-RWJo], 



ClO 
Cll 

C'lo 



R'^rn^ 

R±Vi 
_R!^ 

R± 



xx\x 4- 2x' - 2)Ji, 



x{x -I- x' — 1) Ji , 

2 

-xx\R\^ (1 — x') Jo -f i?^xJi], 



= , 
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^00 



^01 



^00 



-^xx'[R'Aix - 2) Jo + R±{2 - 3a;') Ji], 

m 

R' 

\x\R',{x — 1) Jo — R^_x Ji], 

-{r'\{x - 1)xJq - R'.Rj_xx'Ji, 
h ^ 



- 3) - 3.t' + 4 



Cqi 



;(a;(a;'-3) 

(a; - l){x - l){x'fi^j + a;^^,) Joj, 
-a;(2x + x' - 2)Jo, 



- 1) a;' Jo 



c'n 
Cio 



R' 

-^xx'\R\ xJa - Rdx' - 1) Ji], 

m 



= , 



III I I 

— XX ix + 2x — 2) Jo, 
~-x{x + a;' — 1) Jo- 



Appendix E: Coefficients in the equations (53| 



The three-body Fock component (one fermion + 
two bos ons) is decomposed in four spin structures 
by Eq. (47) with the coefficients 5i_4 being scalar 



functions. These coefficients are finear combinations 



( 53 ) of the functions h and H which are the solution 



of the equations ( 42 ) . The coefficients of these linear 



combinations are given below. 
With the notation 



we have 
an 



m = m^xl + (1 - X2)ii\ + i?2j_, 



m(Ri_L-R2_L)a;2(l + Xi - X2) 
2Kl^xi{l - X2)m 

(Rij_-R2_l) 



2mi?2^(l - 2:2) ' 

?2 



^10 = 

020 = 
^20 = 
^21 = 



■mx2[Rl^{l - X2) + (Ri_l-R2_l)xi] 
2i?2^a;i(l - a;2)?72 

a;2[(l - a:2)(Ri_L-R2_L) + Rl^Xi] 

2mxi{l - X2)'r]2 
mx2{l + xi~ X2) 

2xi{l-X2)'ll2 

1 

2m(l — 2:2) ' 



030 

031 



771^X2(1 + Xi — X2) 

2Rl^xi{\ - X2)m ' 
m 



A 



30 



040 



2i?f^(l-.T2)' 

171^X2 
2i?2^(l-X2)r/2' 

ma;2 
2a;i772 ' 



All = 0-21 = ^31 = 041 = ^40 = Ml = 0. 



Appendix F: Coefficients entering Eq. ( |68[ ) 

The three-body contributions to the form factors 
Ji and F2 are integrals from bi-linear combinations 
of the four spin components 51-4. The coefficients 
determining the three-body contribution to the form 
factor Fi, Eq. (68a), have the form 



^11 



= m4[i?2^ + (Rii-A)(l-a:i)], 



C. 



(1) _ 
22 — 

r<(i) _ 



'-'33 — 
'-'31 — 



-m2{(Ri^.R2i)2 

[(1 - a:i)(R2±-A) - a;2(Ri±-A)](Rii-R2±) 
(Rii-A)i?2^(a:i - 1) + Rl^[{Ii2^- ^)x2 - RI^_]}, 
1 ^(i)r-(i) 

^6^11 W4 . 

m2[(Ri^.A)(Rii.R2i) - i?2^(R2±-A)](a;i - 1). 



The coefficients determining the three-body contri- 
bution to the form factor F2 , Eq. ( 68b I , have the form 



^12 
(^41 

^32 



-4m4(Ri^.A), 
-4m2[(Ri^.A)(Ri^.R2i) 

/^(2) 
L'41 : 



RlA^2±-A)] 



(2) 
34 



r(Rii-A)a 



(1) 
44 • 



We remind that = A^ 
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